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Abstract. We present a simple analytical method to de- 
scribe the structure of a spherically expanding envelope 
with strong mass outflow. The structure is consistently 
connected to the hydrostatic stellar interior and provides 
an adequate description of the outer boundary conditions 
Q«>^ for stellar models with large mass loss rates. 

We apply our treatment to evolutionary models of 
O ^Wolf-Rayet (WR) stars in order to study the possible in- 
fluence of the stellar winds on the interior, and to deter- 
mine more reliable radii of WR stars. Independently of 
\^ the wind parameters (wind density, opacity, velocity law) 
the interior structure and evolution of WR stars is found 
to be unaffected by the outer layers. 
On On the other hand, the stellar parameters (radii, ef- 
C fective temperatures) may well depend on the wind struc- 



OS 



ture. For hydrogen rich WR stars (WNL) we flnd the exis- 



^—^ fence of a temperature domain in the HR-diagram, where 
a transient concentration of stars on their blueward track 
On is predicted in case of a strong backwarming from the 
*^ wind. 

rS I For WNE and WC/WO stars with strong mass loss 
Q rates we also derive subphotospheric radii corresponding 
to Rosseland optical depths of r ~ 10-20. The depen- 
ly^ dence of the subphotospheric radii on the adopted enve- 
lope structure is discussed. With respect to wind-free stel- 
lar models the subphotospheric radii are increased by up 
to a factor of ~ 4 for the most luminous WNE or WC stars. 
These radii and the corresponding effective temperatures 
should roughly be comparable with the stellar parameters 
("core" radii and temperatures) of non-LTE atmosphere 
models of WR stars. Comparisons using the newly derived 
subphotospheric radii yields a better agreement with ob- 
servations. 

The stellar parameters obtained with the new treat- 
ment allow a better assignment of theoretical spectra to 
evolutionary tracks of evolved WR stars (WNE, WC). 
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This also provides the base for future studies of the spec- 
tral evolution of post main-sequence massive stars and 
their descendants. 

We also point out the possible importance of the iron 
opacity peak for the acceleration of WR winds in the op- 
tically thick part, which may be essential for the under- 
standing of the dynamics of WR winds. 

Key words: Stars: evolution - atmospheres - mass-loss 
- Wolf-Rayet - Hertzsprung-Russel diagram 



1. Introduction 

Wolf-Rayet (WR) stars are among the stellar objects 
showing the largest mass loss rates observed (M ^ 
10"'* MQyr""'^). As a consequence several complications 
arise for the derivation of stellar parameters and compar- 
isons with evolutionary models, which mainly limit quan- 
titative comparisons with observations to luminosities and 
surface abundances. The largest difflculty is the determi- 
nation of a stellar radius, which is not unique (cf. Baschek 
et al. 1991). Essentially there are three reasons: i): As the 
photospheres of WR stars can be formed in the supersonic 
part of the wind, the density scale height is large, and 
the radius at a given optical depth becomes wavelength 
dependent (see e.g. Hillier 1987). ii): The density struc- 
ture of the winds is not well known. Although it is widely 
thought that the basic driving mechanism of WR winds is 
radiation pressure (cf. Lucy & Abbott 1993), the dynam- 
ics of these stars is not yet understood (Schmutz 1994, 
1995 and references therein). Therefore, the density struc- 
tures are uncertain and have to be assumed or determined 
empirically, m:) With present state-of-the-art models we 
are still not able to account consistently for all relevant 
opacity sources in the expanding envelope (Schaerer et 
al. 1995c). 



is also of interest in a more general context. Indeed wind- 
free evolutionary models predict an important population 
of WR stars with temperatures ^ 125000 K. If theoret- 
ical spectra from plane parallel atmosphere models cor- 
responding to these temperatures are however assigned 
to this population, very strong nebular excitation is pre- 
dicted, in contrast to e.g. observed H II regions (cf. Garcia- 
Vargas & Diaz 1994). It is, however, also not correct to 
use photospheric temperatures (12/3) to characterise the 
spectral energy distribution (see below). Therefore, to ob- 
tain an adequate description of the emergent UV flux of 
a young stellar population including WNE and WC/WO 
stars requires a good knowledge of the fundamental pa- 
rameters of these stars. This is also fundamental impor- 
tance to determine the role of starbursts in AGN's (cf. Cid 
Fernandes et al. 1992, Leitherer et al. 1992). 

Evolutionary models have so far treated stellar winds 
by, what I here call "a posienon" methods. This means 
that the stellar interior (hydrostatic or hydrodynamic) is 
usually treated with a plane parallel grey atmosphere as 
the outer boundary condition. The spherical extension of 
the atmosphere is neglected, and no backscattering from 
the dense wind (also called wind blanketing) to the inte- 
rior is taken into account - in other words the interior is 
completely unaware of the presence of the wind ! For pho- 
tospheres of OB stars the effect of wind blanketing, which 
leads to backwarming and also modifies the spectral en- 
ergy distribution, has been studied by Abbott & Hummer 
(1985). Standard evolutionary models account for the ex- 
tension of the atmosphere by translating the stellar radius 
and the effective temperature Tgff of the wind-free star, in 
an independent "a posteriori" step, to a photospheric ra- 
dius and corresponding effective temperature (called R2/3, 
and T2/3 respectively), which refer to an optical depth of 
2/3. 

Several approaches to obtain an "a posienon" es- 
timate of i?2/3 and T2/3 have been made, taking ei- 
ther electron scattering opacity into account (de Loore 
et al. 1982, Langer 1989a), or including additional line 
opacity (Schaller et al. 1992), which yields a better agree- 
ment with observations (Maeder & Meynet 1994). On 
the other hand the spectral energy distribution of WR 
stars with identical photospheric temperatures T2/3 can 
differ significantly. It is instead characterised by the so- 
called "core temperatures" corresponding to approx- 
imately unit thermalisation depth, i.e. optical depths of 
10-20 (Schmutz et al. 1992). A unique correspondence be- 
tween a position in the HR-diagram and the stellar spec- 
trum of WR stars thus requires the knowledge of "core" 
radii from evolutionary models. This is one aim of the 
present work. 

The first combined stellar structure and atmosphere 
(CoSiar) models, which consistently treat the stellar in- 
terior and their wind, have recently been presented by 
Schaerer (1995) and Schaerer et al. (1995ab, hereafter Pa- 



also predict the detailed spectral evolution for massive 
stars taking non-LTE effects and line blanketing into ac- 
count. The atmosphere structure for OB stars consists of 
a quasi-hydrostatic spherically extended photosphere and 
the wind. One difficulty, however, arises for WR phases 
(e.g. WNE, WC) with very large wind densities, since sub- 
photospheric layers (r ^ 10) may still be located at su- 
personic expansion velocities. To treat the zone reaching 
further down to the hydrostatic interior (subsonic veloci- 
ties) with their detailed non-LTE radiation transfer calcu- 
lations is clearly prohibitive. For WR stars an alternative 
treatment is thus required. Therefore, we here introduce 
a simple method to determine the temperature structure 
in a spherically expanding envelope, which may be of very 
large optical depth. The density- and velocity structure 
are parametrised and guarantee a smooth transition with 
the hydrostatic interior. Our formulation not only yields 
"core" radii characterising the emergent spectra, but it 
also allows to account consistently for the infiuence of the 
expanding envelope (backwarming) on the stellar interior. 
Due to the simplicity of the approach presented here, it 
can easily be applied to different WR phases, and the small 
parameter space is rapidly explored. 

The outline of the present publication is as follows. 
First we summarise the physical ingredients of our models 
and develop the equations describing the expanding enve- 
lope (Sect. 2). Envelope models are presented in Sect. 3, 
where the infiuence of the modified boundary conditions 
is studied in detail. Consistent evolutionary models for 
WN and WC stars including the mass outfiow are given 
in Sects. 4 to 6, where we also discuss its impact on the 
stellar parameters and the interior evolution. A short com- 
parison with observations is made in Sect. 7. Finally the 
main conclusions are summarised in Sect. 8. 



2. Input physics 

The present calculations are based on the Geneva stel- 
lar evolution code. Taken apart the new treatment of 
the atmosphere and envelope described below, we use the 
same physical ingredients as the latest grids of Meynet et 
al. (1994). For a more detailed discussion of most ingre- 
dients the reader is referred to Schaller et al. (1992). In 
particular the models are computed with a moderate core 
overshooting of d/ Hp = 0.20, where d is the overshooting 
distance, and Hp the pressure scale height at the bound- 
ary of the classical Schwarzschild core. For possible exter- 
nal convection zones we follow the treatment of Maeder 
(1987), which includes turbulent pressure, acoustic fiux, 
and a density scale height (cf. Sect. 3). Radiative opaci- 
ties are from OPAL (Iglesias et al. 1992), completed with 
the atomic and molecular opacities of Kurucz (1991) at 
low temperatures. For WR phases with modified carbon 
and oxygen surface abundances (WC), we use the recent 



Rogers 1993). 

2.1. Spherically outflowing envelope and atmosphere 

In order to consistently model the stellar wind and the 
interior we introduce a simple analytic method to treat 
the outflowing envelope and atmosphere. This treatment 
solves for the mass conservation, and radiative equilibrium 
in a spherically expanding medium. The momentum equa- 
tion, which determines the velocity and density structure, 
is replaced by adopting a parametrised velocity law 
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where t;oo is the terminal velocity. The free parameter /3, 
with typical values of 1-2, expresses the steepness of the 
velocity held. Rq will be determined below. Equation 1 
and the continuity equation 



M = A-Kr^ p(r)v(r) 



(2) 



determine the structure of what is henceforth called the 
expanding envelope. Mostly for convenience we chose the 
more general term "envelope" rather than atmosphere, 
since the total optical depth of this region can be sub- 
stantially larger than what is usually understood as the 
atmosphere. 

In order to allow a smooth transition to the hydrostatic 
interior model, the lower boundary of the expanding en- 
velope is deflned at a given velocity, which, as a simple 
working hypothesis, we choose at a constant fraction s of 
the isothermal sound speed of the perfect gas. Varia- 
tions of s also allow to explore the effect other deflnitions 
of the sound speed (e.g. adiabatic) might have. The inner 
boundary is thus given by 



v(R^) 



where 



(3) 



p is the mean molecular weight, and R^ and are the 
radius and effective temperature of the hydrostatic star re- 
spectively. From the continuity equation (Eq. 2) this con- 
dition determines the density at the inner boundary, as 

well as Ro = R^ |l — oo)^''''| required in Eq. 1. 

The temperature T(R^) is calculated from the analyt- 
ical solution of the spherically extended grey atmosphere 
in radiative equilibrium using the generalised Eddington 
approximation derived by Lucy (1971; cf. also Lucy & Ab- 
bott 1993). The temperature structure is given by 
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where W is the geometrical dilution factor deflned by some 
photospheric radius Rp > R^. Note that the optical depth 



low). To clarify the non-unique deflnition of an "effective 
temperature" Tgff in objects with spherical extension, we 
repeat that it is deflned from the total emergent luminos- 
ity L = Airr^aT^ff. Thus, the effective temperature of 
the hydrostatic star is e.g. given by setting r = R^ and 

= Tgff. In the same way the effective temperature can 
be deflned at any depth r. In particular T2/3, Tio, . . . de- 
note henceforth the effective temperature at given optical 
depths T = 2/3, 10 . . . 

With the adopted velocity law, the effective optical 
depth f = J^°° Kp (■^) dr can be integrated analytically 
for a depth-independent opacity. For values of /3 7^ 1, 2, 
and 3 this yields 



t(x) = C [{p - ■ 
(l_^)-/5+2] 
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where 
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and X = Ro/r. Note, however, that the optical depth r = 
Kpdr is given by 

r(a;) = C(l-/3)-i [l-(l-a;)^-''] for /3 7^ 1. (7) 

Unless specifled otherwise, the opacity k is taken as the 
electron scattering opacity ag of a fully ionised mixture 
with abundances from the stellar interior model at the 
given stage. 

Since we are dealing with optically thick envelopes one 
has W = 1/2 at the inner boundary R^, and the tempera- 
ture T(R^) is thus solely determined by f. The variable C 
describes the behaviour of f on the basic wind parameters 
(M, Voo, opacity, and radius R^). From this we immedi- 
ately see that for a given velocity law the temperature at 
the inner boundary R^ increases with mass loss, increas- 
ing opacity, decreasing terminal velocity, and decreasing 
stellar radius, as expected. 

The dependence of the temperature on the shape of the 
velocity law is illustrated in Fig. 1 for two different values 
of v^/voo- The solid line (v^/voo = 0.01) represents a typ- 
ical value for WN stars, if we set equal to the sound 
speed a^. For comparison, a case where the inner bound- 
ary is deflned at lower velocity is also shown (dashed). 
Figure 1 shows that for all values of /3 ^ 1 the depth scale 
t(R^), and hence the temperature at the inner boundary 
of the wind does not strongly depend on the shape of the 
velocity law. For v^/voo = 0.01 e.g., the temperature at 
the bottom of the wind varies by less than ~ 15 %. The 
changes resulting from the adopted value of v^, and the 
uncertainties in the opacity k are clearly larger than the 
exponent of the velocity law (cf. Sect. 3). 
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Fig. 1. Logarithm of f/C at the inner boundary i?^ as a func- 
tion of the exponent /? of the velocity law (Eq. 1). Ro is implic- 
itly defined by the choice of Vi^/vco. The solid line correspond 
to the typical value of v^^ = for WN stars, where ~ 20 
km s~^, and v^d ~ 1000-2000 km s~^. Note that for all values 
of /? > 1 the depth T{Ri^) varies by less than ~ 0.23 dex, which 
corresponds to temperature changes of less than ~ 15 % 



Once the temperature T(R^) is determined (p^ is also 
known), we finally compute the total pressure at R^, which 
is also required for the integration of the stellar structure 
equations. The pressure is derived from the equation of 
state taking into account partial ionisation of H, He, C, O, 
Ne, and Mg (cf. Schaller et al. 1992). It is straightforward 
to replace the usual boundary conditions of stellar interior 
models (plane parallel grey atmosphere, cf. Kippenhahn 
& Weigert 1990) with the method presented above. The 
only additional parameters, which need to be specified, are 
the wind velocity law (t;oo , /3) and the wind opacity. The 
terminal velocity t;oo can be obtained from observations. 
The mass loss rate has already been taken into account 
to determine the mass removed at each evolutionary time 
step. 

With this procedure we have thus obtained an ade- 
quate method to describe the outfiowing envelope and at- 
mosphere with few basic parameters, which describe the 
hydrodynamic structure from close to sonic velocities up 
to the highly supersonic asymptotic wind velocity. Assum- 
ing the velocity structure is monotonous in the entire star, 
the choice of the velocity at the inner boundary guaran- 
tees that below hydrostatic equilibrium should be closely 
verified, which means that hydrostatic interior solutions 
are clearly appropriate. 



Before turning to the remaining physical ingredients of 
our models, it is worth clarifying some problems relating 
to the radius and temperature definition of the extended 
atmospheres of WR stars (see also Hillier 1991, Schmutz et 
al. 1992, Schmutz 1995). As already mentioned the major 
difficulties stem from the hydrodynamic structure (espe- 
cially at large optical depths) and the wind opacity, which 
are badly known. For the first point we have adopted 
a parametrised description as a working hypothesis (see 
Sect. 2.1). A large uncertainty in our treatment arises from 
the difficulty of specifying the wind opacity. In this context 
it is, however, important to clearly distinguish the two fol- 
lowing points: (i) the cooling and heating processes, which 
determine the temperature structure (including the back- 
warming), and (ii) the (frequency dependent) continuum 
opacity allowing to define the radius, where radiation is 
thermalised. Let us therefore discuss these points sepa- 
rately. 

(i) As will be shown in Sect. 4, the effect of backwarming 
can have an infiuence on the outer layers of WN stars 
with hydrogen, and lead to a relative increase of the 
stellar radius. The large opacity adopted in Sect. 4 (20 
times the electron scattering opacity), corresponds to 
the typical value proposed by Lucy & Abbott (1993), 
which is derived from the effective scattering coefficient 
((Teff) required to explain the dynamical wind struc- 
ture by radiation pressure. (Teff is thus a fiux weighted 
mean opacity, which shows that the Lucy & Abbott ap- 
proach for the determination of the temperature struc- 
ture is basically equivalent to the treatment of Schaller 
et al. (1992) to include the effect of line blanketing. 
Whether, however, the procedure using fiux weighted 
mean opacities indeed yields a good description of 
the detailed heating and cooling mechanisms includ- 
ing line blanketing is not clear. The present method of 
Schaerer & Schmutz (1994) applied to WR stars gives 
only a modest temperature increase (cf. Schaerer et 
al. 1995c). This question awaits a quantitative treat- 
ment including the calculation of thermal and statisti- 
cal equilibrium. 

(ii) Given a strong variation of the continuum opac- 
ity with wavelength, any radius definition in ex- 
tended atmospheres is necessarily wavelength depen- 
dent (e.g. illustration in Hillier 1987, Fig. 2). For WNE 
stars e.g., the electron scattering depth Tg = 2/3 yields 
a good estimate of the thermalisation optical depth 
Tth = 2/3 at UV to visible wavelengths (see Schmutz 
1995), and hence for the photospheric radius character- 
ising the largest fraction of the directly observable fiux. 
On the other hand, to obtain a radius completely in- 
dependent of frequency requires radiation at all wave- 
lengths to be thermalised, which is typically the case at 
optical depths tross ~ 10-20. Referring to these opti- 
cal depths therefore allows the physically most relevant 



trum including the EUV flux shortward of the Lyman 
edge. 

For the non-LTE conditions in WR winds the Rosse- 
land opacity is mostly dominated by electron scatter- 
ing. Therefore the optical depths, and hence the radii 
R(t ^ 10) predicted in this work, should be similar to 
those from non-LTE model atmospheres, provided the 
other parameters are identical. Changes of the opac- 
ity due to the ionisation stratification, are not taken 
into account in our procedure, but should be of minor 
importance. 

To complete the physical ingredients of our models we 
shall now turn to the description of the wind parameters. 

2.3. Mass loss rates 

Since mass loss is of particular importance for WR phases, 
the adopted prescriptions shall briefiy be recalled here. 
Remember that in our evolutionary models the entry in 
the WR phase is defined by a hydrogen surface abundance 
X < 0.4 in mass fraction, and Tgff > 10000 K. 

— For the WNL phase we take an average mass loss rate 
of 81O-5M0yr-i, as in Meynet et al. (1994). Note 
that the observational scatter is found to be ~ 0.4 
dex, with a mean being slightly lower than the adopted 
value (Hamann 1994, Crowther et al. 1995b). 

— For WNE and WC/WO stars, we use the mass de- 
pendent mass loss rates from Langer (1989b), as in 
previous evolutionary grid calculations (Schaller et 
al. 1992, Meynet et al. 1994). The adopted mass loss 
rates show a good agreement with the observations 
presently available (cf. Hamann 1994). 

2.4. Wmd velocities 

In this paragraph we describe the choice of the adopted 
terminal velocities for WR stars. 

2.4.1. WN stars 

Observed terminal velocities of WN stars are in the range 
of ~ 700 to 2500 km s"-'^. Earlier subtypes tend to have 
higher t;oo although the relation shows a large scatter 
(e.g. Hamann et al. 1993, Crowther et al. 1995b). For WR 
stars no clear picture has emerged yet, which would re- 
late the terminal velocity to fundamental stellar param- 
eters. In analogy with OB stars, one can suspect a cor- 
relation of Voo with the escape velocity. Using the stel- 
lar parameters of WNL stars from evolutionary models 
(cf. Meynet et al. 1994) and applying the same propor- 
tionality as observed for O stars (t;oo = 2.78i'esc, Groe- 
newegen et al. 1989), this yields in fact values, which are 
compatible with the observations. On the other hand, one 
is obviously tempted to relate the trend of higher t;oo for 



crease of t;oo towards larger temperatures, which is found 
in the data of Hamann et al. (1993), can be fitted with 
the following formula: 

= (2395.87 ± 244.83) logT^ - (9656.98 ± 1125.04), (8) 

where the residual rms is 241.61 km s~^. At temper- 
atures logT^ < 4.24 we adopt t;oo = 500 km s""'^, 
which corresponds to typical values for Ofpe/WN9 stars 
(e.g. Crowther et al. 1995a). For the purpose of this work, 
Eq. 8 is clearly sufficient to reproduce the observed trend. 
For most calculations we shall however simply use a con- 
stant value of Vco = 1400 km s""'^, which is close to the 
mean value of 1347.5 km s~^from the data of Hamann et 
al. (1993). As will become clear later (see Sect. 4), the ex- 
act choice of the terminal velocity in the WN phases is, 
however, not crucial for our results. 

2.4.2. WC stars 

For Wolf-Rayet stars from the WC/WO sequence, a clear 
correlation of the terminal velocity with spectral class is 
well established (cf. Abbott & Conti 1987) and has been 
been confirmed by recent studies, which take turbulence 
into account to obtain accurate line profile fits (Eenens 
& Williams 1994). On the other hand there is also am- 
ple evidence for a correlation of the C/He abundance ra- 
tio with WC/WO subtype (Smith & Hummer 1988, Ee- 
nens & Williams 1992) and possibly also for 0/He (Nugis 
1991). Adopting the relation between (C-|-0)/He and the 
WC/WO subtype from Smith & Maeder (1991), one is 
thus able to assing a wind velocity provided the abun- 
dances of He, C, and O are known. 

The method is the following: We adopt the terminal 
velocities from Eenens & Williams (1994) derived from 
IR lines taking turbulence into account. The observations 
include 18 WC stars, and cover the largest subtype range 
for which homogeneous data is available. Their values are 
plotted in Fig. 2 (triangles). They can well be fitted by 

= (3738.49 ± 236.18) - (281.40 ± 31.07)WC, (9) 

with a residual rms of 190.85 km s~^. WC stands for the 
subtype number, and t;oo is to be taken in km s~^. As a 
comparison we have also plotted the values derived from 
UV lines by Prinja et al. (1990, squares) together with 
a fit to their data, which show a slightly steeper trend. If 
the latter relation is extrapolated to WC4, we see that the 
agreement between both fit relations is of the order of ^ 10 
% for Voo, which is clearly sufficient for our purposes. The 
terminal velocities for WC stars with strong lines (WC- 
s, Koesterke & Hamann 1994) agree well with the above 
relation, while their weak lined WC stars (WC-w) have 
considerably larger wind velocities. 

From Smith & Maeder (1991) we now take the mean 
values of (C-|-0)/He (in number) to be representative for 
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Fig. 2. Observed terminal velocities for WC stars as a function 
of subtype. Triangles: 18 stars from Eenens & Williams (1994), 
open squares: 16 measurements from Prinja et al. (1990). The 
solid line shows the fit of Eq. 9 to the values of the former 
authors, the dotted line a fit to the Prinja et al. data 



region treated here covers only a very small percentage of 
the total mass, its extension corresponds to ~ 60 - 90 % 
of the outer radius. We solve for the continuity equation, 
hydrostatic equilibrium, and the energy transport allow- 
ing both for radiative or convective transfer. In a next step 
(Sects. 4 and 6) we will then present evolutionary tracks, 
where the entire stellar structure is solved consistently. 

3.1. Calculated models 

To illustrate the envelope structure and its dependence 
on different boundary conditions we have chosen a model 
corresponding to the WNL phase from the 85 Mq track 
calculated with our "combined stellar structure and at- 
mosphere models" (Schaerer 1995, Schaerer et al. 1995a). 
The stellar parameters are summarised in Table 1. The 
temperature corresponds to a domain, where WNL stars 
are found to be sensitive to boundary conditions ("WNL 
strip" see Sect. 4). Qualitatively the results do not de- 
pend on the adopted luminosity. Therefore the discussion 
in this paragraph also applies to WN models of lower ini- 
tial masses located in the "WNL strip" . 



Table 1. Adopted stellar parameters for envelope models 



the corresponding WC subclass, i.e. WC9: 0.045, WC8: 
0.1, WC7: 0.2, WC6: 0.3, and WC5: 0.55. A linear least 
square fit yields (C + 0)/Re = 1.086 - 0.121 WC, where 
again WC is the subtype number. Combined with Eq. 9 
we can express the terminal velocity as a function of 
(C-|-0)/He, which yields (in kms~^) 

'"'^''^ — C + O\(10) 



He 



3738.49-2324.36 1.086 



He 



This formula will be used to determine t;oo in WC/WO 
phases. Despite its simple derivation it reproduces well 
the observed trend of increasing wind velocities for WC 
stars, and also yields a good agreement with observed t;oo 
values of WO stars (Kingsburgh & Barlow 1991), when 
Eq. 10 is extrapolated towards larger (C-|-0)/He values 
than formally allowed. 

3. Envelope models the influence of the mass out- 
flow 

In this paragraph we construct envelope models including 
the spherically expanding stellar wind and the lower hy- 
drostatic layers. This allows us to study the infiuence of 
the backwarming effect of the wind on the external lay- 
ers of the stellar structure. An envelope model describes 
the outermost stellar layers and is defined by the outer 
boundary conditions (cf. Table 1), the atmospheric bound- 
ary conditions (Sect. 2.1), and the mass it embraces (here 
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The infiuence of the atmospheric boundary condition 
on the envelope solutions can be illustrated with several 
models, which are summarised in Table 2. The first model 
has been computed with a plane parallel grey atmosphere 
in the Eddington approximation. The inner boundary of 
the atmosphere is set at tross = 2/3. The remaining mod- 
els (2-4) have been calculated with the outfiowing atmo- 
sphere models from Sect. 2.1. Their inner boundary is de- 
fined by the velocity t;^, for which we take values no larger 
than the isothermal sound speed to guarantee a smooth 
transition with the hydrostatic part. Models 2 and 3 ex- 
plore the infiuence of the opacity in the wind. The value 
for model 2 corresponds to the electron scattering (Tg for 
the entirely ionised mixture. In model 3 we enhance the 
scattering opacity to simulate strong line blanketing. The 
adopted enhancement factor of 20 is consistent with the 
effective scattering coefficient required if radiative driving 
dominates the dynamics of WR winds (cf. Lucy & Abbott 
1993). Finally, in model 4 the inner boundary is modi- 



Table 2. Summary ol envelope models, r mdicates plane parallel atmosphere, W wmd. ihe mput parameters are p, and 
the wind opacity k (cf. 2). The other quantities denote the value of the optical depth, temperature, pressure, and density at the 
inner boundary of the atmosphere 
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Fig. 3. Density structure of the envelope models 1 to 4, plotted 
as a function of log P used as the depth variable. The symbols 
are as follows: model 1 (short-dashed), 2 (solid), 3 (dotted), 
and 4 (long-dashed). Units on all structure plots (Figs. 3 to 7) 
are in cgs 



3.2. Discussion 

In Figures 3 and 4 we have plotted the run of the density 
and radius of all models. The total pressure serves as a 
depth variable^. The model with the plane parallel atmo- 
sphere providing the outer boundary conditions (model 

^ Note that since radiation pressure strongly dominates, and 
the temperature gradient is essentially independent of P and 
T in the entire envelope, all temperature profiles would lie on 
the same line if plotted as a function of the pressure. They are 
therefore not shown here. 
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Fig. 4. Radial extension of the envelope models (same symbols 
as in Fig. 3). All models have a given external radius (cf. Table 
1). The inner boundary is set at Mr/M = 0.997. Note in par- 
ticular the difference of the r gradient in outer parts of the 
envelope, which leads to the different inward extensions 



1), and the wind model 2 show essentially the same struc- 
ture. Despite the small temperature increase in model 2 
(cf. Table 2) they are nearly indistinguishable in these fig- 
ures. For the parameters of model 2 the wind effects are 
therefore negligible. 

Both models 3 and 4 have a higher pressure at the 
outer boundary due to the strong "backwarming" from 
the wind. The case differing most from the plane parallel 
atmosphere is clearly model 4, since it also has a larger 
density. How does this infiuence the depth of the enve- 
lope ? From Fig. 4 we see that the models show different 
inward radial extensions. This behaviour is simply under- 



be written as 
dlogr rP 

dlogP ~ ~ GMrP ' 

in Lagrangean coordinates. In the envelope is nearly 
constant, and the outer value for the radius is fixed (r = 
Ri,). Therefore both a higher pressure and lower density 
imply a steeper r gradient, and lead to a larger exten- 
sion of the envelope. Considering variations of the outer 
boundary conditions, a large envelope extension can thus 
be obtained by increasing the ratio of P / p. If the density 
at the outer boundary can be kept constant, a tempera- 
ture increase (due to larger opacity, or a change of [3 as 
illustrated in Fig. 1) at the bottom of the atmosphere thus 
leads to a larger inward extension. This explains the differ- 
ences between models l-\-2 and model 3, which are, how- 
ever, rather small. The main reason for the small radius 
difference is that in model 3 the density rapidly converges 
to the solution of the "classical" model 1. 




log P 

Fig. 5. Detailed view of the outer zones of models 1 (plane 
parallel atmosphere) and 4 (wind). The shaded area denotes 
the convective zone of model 1, which is due to the Fe opacity 
peak. In model 4 the convection zone starts at its outer bound- 
ary and has the same inward extension as model 1. Shown is the 
density structure as a function of the total pressure. Structure 
of model 1 (solid), and model 4 (long dashed) calculated with 
a density scale height and II Hp = 1. Calculations adopting 
the pressure scale height yield density inversions: dotted line 
(model 1), dashed-dotted line (model 4). The short-dashed line 
shows the temperature (scale on the right) 
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Fig. 6. Detailed view of the outer zones of models 1 and 4. 
The shaded area is as in Fig. 5. Shown are the Rosseland mean 
opacity (in cm^ g~^)) ^nd the total local luminosity Lr in Ed- 
dington units for model 1 (solid), and 4 (long dashed). The 
short-dashed line shows that the radiative luminosity of model 
1, using the density scale, never exceeds the Eddington lumi- 
nosity. Discussion in text 

To understand the important increase of the extension 
of model 4 it is worth recognising an additional "fortu- 
itous" complication, which turns out to be relevant for 
determining stellar radii of hydrogen rich WR stars. We 
have therefore plotted an enhanced view of the outer en- 
velope structure of the models with the largest differences 
(models 1 and 4) in Figs. 5 and 6. The point of particu- 
lar importance is the finding of a small convective zone, 
which is present in all our models. The convection zone, 
also noted by Glatzel & Kiriakidis (1993), Alberts (1994), 
and Schaerer et al. (1995c), is due to the opacity peak of 
iron (cf. Iglesias et al. 1992), situated around logT ~ 5.2 
(see Fig. 6). 

The key point to explain different extensions is that 
in the convective zone the density increases less steeply 
than in the radiative layers. The reason for this will be 
discussed below. If we consider e.g. model 4, we see that 
due to the convective layer, a lower density is maintained 
between 6 ^ log_P ^ 6.8 compared to models 1-3. There- 
fore, according to Eq. 11, the radius of model 4 decreases 
more rapidly in this interval, which leads to the final ~ 40 
% difference in radius. 

Let us now examine the behaviour of the density struc- 
ture within the convection zone. The approximately con- 
stant slope in the \ogp-\ogP diagram (Fig. 5) is due 
to adopted treatment using a density scale height and 
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Fig. 7. Run of the radius as a function of the total pressure de- 
pending on different treatment of convection. The shaded area 
denotes the convective zone of model 1 (cf. Fig. 5). As reference 
models we plot model 1 (solid) and 4 (dashed) both computed 
with II Hp = 1. The long-dashed-short-dashed line is model 
1 but adopting l/Hp = 0.5. The models showing the largest 
extension are those showing the density inversion (I/Hp = 1.6 
model 1 and 4, dotted and dashed-dotted lines, nearly coincid- 
ing) 

l/Hp = 1 as in Maeder (1987). If instead we use the pres- 
sure scale height and l/Hp = 1.6, we obtain a density 
inversion, which occurs where the radiative acceleration 
exceeds the local gravity g, i.e. K_Lrad/(47rr^c) > g. This is 
illustrated in Fig. 6 for models 1 and 4. We note that the 
finding of a density inversion is not due to a hydrostatic 
treatment. As mentioned above, for MS stars Glatzel & 
Kiriakidis (1993), using a hydrodynamic code, also find a 
/9-inversion related to the Fe opacity peak (cf. Gautschi 
& Glatzel 1990). This may, however, also be related to 
the use of time steps, which are larger than the relevant 
timescales (Meynet, 1994). 

Before we discuss its impact on the extension of the 
envelope, it is worth clarifying the question of density in- 
versions in this context (see Maeder 1992 for a recent re- 
view). First, we note that the region with the p-inversion 
is not Rayleigh- Taylor unstable, since there the effective 
gravity is directed outwards (cf. Wentzel 1970). Similarly, 
the strong supra-Eddington luminosity observed within 
the convection zone (Fig. 6) is simply due to the fact that 
other modes of energy transport contribute with a sig- 
nificant fraction to the total energy transfer (up to ~ 20 
% convective, and 40 % acoustic fiux respectively in the 
considered cases). If the density inversion is suppressed 



minosity obviously never exceeds the Eddington limit. In 
other words, the existence of the convective zone is due to 
the strong increase of the opacity, which, according to the 
treatment may or may not lead to a density inversion. The 
(total) supra-Eddington luminosity is related to the strong 
contribution of non-radiative fiuxes, and does, however, 
not necessarily lead to an outward acceleration of mat- 
ter, as suggested in particular by Kato & Iben (1992) for 
WR stars. Nevertheless, we note that in the outer layers 
the effective gravity is indeed strongly reduced, indicating 
that "potentially" there is a large radiative force avail- 
able, which might contribute to initiate the mass outfiow 
(cf. Sect.4). 

Given the dependence of the total envelope extension 
on its outer density structure, and hence on the treatment 
of convection, we have also investigated its infiuence on 
models 1 and 4. In Figure 7 we have plotted the run of the 
radius using either the density scale height (and l/Hp = 1 
or 0.5) or the pressure scale with l/Hp = 1.6. As expected, 
the models having the lowest density in the outer parts 
show the largest radial extension. Using the pressure scale 
height the wind model with strong backwarming (model 4) 
in fact turns out to be nearly identical to the plane parallel 
case. For the present model parameters, changes due to 
the treatment of convection and different atmospheres are 
thus of similar order. 

Let me summarise. We have seen that for typical ef- 
fective temperatures of WNL stars envelope solutions are 
sensitive to the treatment and changes of atmospheric 
boundary conditions. Changes are indeed expected be- 
tween plane parallel atmospheres, which are certainly not 
adequate for stars with strong stellar winds, and a more 
consistent treatment including an expanding atmosphere 
with wind blanketing. The basic result is that the spherical 
extension of the envelopes may be increased. This finding 
is intimately related to the presence of a small layer with a 
considerable opacity increase (due to iron), which is unsta- 
ble to convection. Depending on the effective temperature, 
i.e. the depth where this layer is located, its importance for 
determining the envelope extension and hence the stellar 
radius is expected to vary crucially. This will be studied 
in the following section, where we present evolutionary 
models covering the WN phase. 



4. WN stars with expanding atmospheres 

Let us now have a look at evolutionary tracks including 
spherically expanding atmospheres. We shall first discuss 
WNL stars, which are defined by the presence of both 
hydrogen and helium on their surface. The more evolved 
WNE and WC/WO phases are covered in Sects. 5 and 6. 

To study the infiuence of the atmospheric treatment we 
have calculated two evolutionary tracks for initial masses 



Note that all temperatures refer to the effective tempera- 
ture of the hydrostatic star (cf. Sect. 2.1). For the WNL 
phase considered in this paragraph, however, the extension 
of the layers between the hydrostatic interior and elec- 
tron scattering optical depths of r = 20 or 10 is small. 
Therefore we have T20 ~ Tiq. The largest difference 

amounts to ~ 0.1 dex (see below). 
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Fig. 8. HR diagram for WR phases of the 40 and 85 
Mq model. Remember evolution proceeds bluewards and to 
lower luminosities. Solid lines show tracks with plane par- 
allel atmospheres, dashed lines with spherically expanding 
envelopes/atmospheres. For strong wind blanketing one ob- 
tains the dotted tracks. Shaded areas indicate the position 
of the WNE and WC phases at log ~ 5.1, and the sug- 
gested strip for WNL stars predicted with strong blanketing 
(marked "WNL"). The ZAMS from Meynet et al. is also plot- 
ted (long-dashed) 

In Figure 8 we have illustrated several tracks for 
each initial mass. First, the dashed line shows the re- 
sults from standard models using the plane parallel atmo- 
sphere. Adopting the spherically expanding atmosphere 
with = a^, l3 = 2.1, t;oo from Sect. 2.4, and k = ag, 
yields the solid line track. We immediately see that the 
differences are negligible for the WNL phases. This is ba- 
sically due to the "moderate" wind density and low opac- 
ity, resulting in relatively small total optical depths. The 

^ The pre-WR tracks correspond to Meynet et al. (1994) for 
the 40 Mq model; the 85 Mq track is from the "combined 
stellar structure and atmosphere models" of Schaerer (1995) 
and Schaerer et al. (1995a). 
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Fig. 9. Time evolution of the temperature during the entire 
WNL phase of the initial 40 Mq model. Tracks are plotted with 
the same symbols as in Fig. 8. The hydrogen surface abundance 
corresponding to both models is shown by the long-dashed and 
dotted lines (scale on the right) 



largest values obtained are t(R^) ~ 10-30 for logT^ = 5., 
which implies a low backwarming from the extended at- 
mosphere. This is in line with the results of the envelope 
models 1 and 2 (Sect. 3). 

On the other hand, if we adopt a larger wind opacity 
(k = 20 (Tg), the backwarming effect is considerably en- 
hanced (cf. models 2, 3). The resulting evolutionary path 
is shown by the dotted line in Fig. 8. Overall the evolution- 
ary paths, and in particular also the deep interior and its 
properties (lifetimes, internal abundance profiles etc.), re- 
main unaffected. Nevertheless one very interesting feature 
can be observed from the tracks with strong wind blanket- 
ing. During the bluewards evolution in the WNL phase the 
models are temporarily stopped, or concentrated in rather 
small strip at temperatures 4.4 <, logT^ ^ 4.65, which has 
been marked in Fig. 8. The timely behaviour is illustrated 
in Fig. 9. While most models advance continuously to- 
wards the blue, the model with strong backwarming main- 
tains a larger radius, thereby staying at ^ 44000 K, 
before it transits in a short lapse of time to larger temper- 
atures. The same behaviour is also found for the 85 Mq 
track, and is luminosity independent. It rests, however, 
on an strong backwarming effect, which is obtained by 
adopting a large opacity in our calculations. As expected 
(cf. Sect. 3), modifying the other parameters is found to 
be of less infiuence. 



from the considerations of Sect. 3. In this temperature 
domain, the small convective zone related to the Fe opac- 
ity is sufficiently close to the surface so that the pressure 
and density is still affected by the external boundary con- 
ditions. In this case the buffering high opacity zone is ca- 
pable to maintain a large radius, as explained in Sect. 3. 
The blueward trend, resulting from the homogeneisation 
due to mass loss, can however not be stopped. Once the Fe 
convection zone has moved sufficiently outwards, density 
and pressure profiles rapidly converge to the same struc- 
ture, which is nearly independent of the outer boundary 
conditions. Thus the bluewards transition occurs. 

As shown in Sect. 3, the importance of the convection 
zone for the envelope extension also depends on the treat- 
ment of the possible density inversion. Calculations using 
either the density- or the pressure scale height, however, 
yield the same "strip" , and show the same dependence on 
the other parameters as discussed above. A brief compar- 
ison with observations of WNL stars in given in Sect. 7.1. 

4.1. Radiation driven winds of WR stars - the role of Fe 

We have seen that the strong opacity peak at logT ~ 5.2, 
due to iron, may play an important role in determining 
the structure of the outermost layers and thereby the stel- 
lar parameters of WNL stars (cf. Sect. 3). In all cases the 
Fe opacity leads to a strong decrease of the local effec- 
tive gravity. In our WNL models this domain is located in 
subsonic layers, whereas along the entire WNE and WC 
sequence (see below) the zone with the opacity peak is 
encountered in the expanding envelope, i.e. at t; > t;^. 
By specifying the envelope structure and adopting a con- 
stant opacity it has therefore completely been neglected. 
As will be shown below this has no implications on the in- 
terior structure and evolution. However, we strongly sus- 
pect that the large Fe opacity bump may be an impor- 
tant contributor to the acceleration of the stellar winds 
of hot WR stars. Similarly it may also be the most im- 
portant blanketing agent in the deep photospheric layers 
(cf. Schaerer et al. 1995c). Observationally Fe V and VI is 
known to be responsible for strong blanketing effects in 
WNE stars (cf. Koenigsberger 1990), while higher ionisa- 
tion stages will dominate at larger depths. A correct treat- 
ment of Fe in line blanketed atmosphere models is however 
not yet possible (Schaerer et al. 1995c). More quantitative 
work will be required to understand the optically thick 
winds of WR stars. Implications of the above suggestion 
(metallicity dependence !?) and the search for possible ob- 
servational constraints also deserve futures studies. 

5. WNE stars 

Let us now study more evolved WR phases. In this para- 
graph and the following ones, we have adopted the fol- 
lowing parameters to describe the expanding envelope: 



of these parameters will be discussed. First we shall con- 
centrate on the radii of WNE stars. The HR diagram from 
the 40 and 85 Mq tracks for the advanced WR phases is 
presented later (Fig. ??). 

After the evolution through the hydrogen rich WNL 
phase, evolutionary models predict a phase showing an es- 
sentially pure He surface abundance. In this phase, identi- 
fied with WNE, and the subsequent WC/WO phases, the 
basic parameters of WR stars follow several simple rela- 
tions (cf. Schaerer & Maeder 1992). The fundamental link 
is the well known mass-luminosity relation (Vanbeveren 
k Packet 1979, Maeder 1983, Langer 1989a, Schaerer k 
Maeder 1992). In particular, a well defined mass-radius 
relation is obtained, independently of whether the stellar 
interior is assumed to be in hydrostatic equilibrium or a 
standard hydrodynamic treatment is accounted for (see 
Fig. 10). 
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Fig. 10. Mass-radius relation for WNE stars. The solid line 
shows the radius i?^ of the hydrostatic interior. The subphoto- 
spheric radii Rio and R20, adopting /? = 2.1, are shown by the 
short-dashed and dotted line respectively. At log M ~ 1.5 the 
same symbols illustrate the radii Rio and R20 obtained with 
/? = 1.1. The thin lines show the relations from a hydrody- 
namic model (dashed-dotted; Langer 1989a) and hydrostatic 
interior models (long dashed; Schaerer & Maeder 1992) 

Sofar all models have used plane parallel atmospheres 
for the outer boundary, which does not allow to take the 
strong mass outfiow into account. Given our treatment 
of the outfiowing envelope we can now also derive its ra- 
dial extension. Since commonly radii of WR stars are de- 
fined at a Rosseland optical depth of typically 10 to 20 



LTE atmospheric models, we have plotted the radius cor- 
responding to these depths (Fig. 10). However, as a word 
of caution we remind the reader that direct comparisons 
with radii derived from non-LTE analysis may still be 
hampered if the hydrodynamic structures (particularly at 
large depths) adopted in both models differ significantly. 
Note that the velocity at r = 20 (10) is between 30 and 
470 km s~^ (100 and 920 km s~^) for the interval plotted 
in Fig. 10. The hydrostatic layers are reached at optical 
depths of t(R^) ~ 20 to 130. Figure 10 shows that for 
WNE stars the subphotospheric layers (i.e. r > 10 or 20) 
may reach up to ~ 2 times the radius of the hydrostatic 
part, which corresponds to a temperature decrease of ^ 
0.15 dex, with respect to temperatures of logT^ ~ 5.1 
- 5.15. The predicted decrease of the spherical extension 
with mass is simply due to the adopted mass-dependent 
mass loss rate, which implies a decrease of the wind den- 
sity with decreasing mass, i.e. luminosity. 

The definition of i?^ obviously depends on the envelope 
parameters (/3, t;oo , k), but also on M, which is however 
not treated as a free parameters (Sect. 2.3). As mentioned 
earlier, we note that the observations show a good agree- 
ment with the adopted mass loss rates (Hamann 1994), or 
even point to a somewhat higher values for M. As we can 
see (Fig. 10) a lower /3 yields slightly smaller radii. The 
terminal velocities, according to Eq. 2.3, are of the order 
of 2550-2680 km s"-'^, which is ~ 15 % larger than typical 
observed values of WNE stars. Most important is clearly 
the wind opacity, for which we have taken the electron 
scattering opacity only. This should represent the mini- 
mum value for frequency independent continuum opacity. 
Provided the hydrodynamic structure of WR envelopes 
does not differ significantly from our description more re- 
alistic radii for WNE stars are thus quite probably under- 
estimated by the mass-radius relation presented here. 

The properties of the hydrostatic interior deserve one 
additional comment. Therefore we have also plotted the 
radius of the hydrostatic interior. As expected this 
agrees well with previous calculations, although we now 
take the backwarming effect of the expanding envelope 
into account. This means that contrarily to WNL stars at 

<, 44000 K, for WNE stars (the same holds for WC/WO 
stars, see Sect. 6) the outer boundary conditions have no 
influence on the interior structure, and in particular on its 
extension. This result is robust to any envelope parameter 
change, since all solutions rapidly converge in the radiative 
layers - the region with opacity peak is now located in the 
expanding envelope (cf. Sect. 4). 

6. WC stars 

In continuation from WNE stars discussed above, we now 
turn to WC stars. Both the wind velocities of WC stars, 
and their mass loss rates differ from the WNE phases. 
While the latter are a factor 5/3 larger for the same mass 



Eq. 9 are smaller during most of the WC/WO phase than 
their WNE counterpart. Thus the winds are denser and 
the extension is increased. 
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Fig. 11. Mass-radius relation for WC/WO stars including the 
spherical mass outflow. Same symbols as in Fig. 10. For each 
radius (Rio, R20) a pair of curves illustrates a variation of 
P from 2.1 (larger radius) to 1.1 (smaller R). Also shown is 
the relation for hydrostatic wind-free stars from Schaerer & 
Maeder (1992; thin long-dashed line). The dashed-dotted line 
illustrates the small radius (i?*) increase due to CO opacities 



This is illustrated in Fig. 11, where we show the mass- 
radius relation for WC/WO stars taking into account the 
spherically expanding envelope. The behaviour is as for 
WNE stars. We note that even with the low wind blan- 
keting adopted, the radii corresponding to r = 10 are 
enlarged by up to a factor of 4 for the models with the 
largest mass loss rate. For this case the total optical depth 
of the expanding envelope amounts to r ~ 160-250 for /3 = 
1.1-2.1. The main effect causing the decrease of spherical 
extension is due to the decrease of M cx M^^^, while t;oo 
only increases by ~ 0.4 dex over the interval shown in 
Fig. 11. Also shown is the small infiuence of the parame- 
ter l3. As for the WNE stars, we would like to point out 
that the radii obtained from this relation should provide 
a lower Urrnt to the "core radii" derived from non-LTE 
analysis, since these are defined on the Rosseland optical 
depth scale by tross ~ 10-20, while our models only ac- 
count for electron scattering opacity. The same words of 
caution (cf. Sect. 5) regarding the hydrodynamic structure 
of both non-LTE models and our approach are, however. 



1 , I 1 I I , , r 




5 4.5 4 

log T 

Fig. 12. HR-diagram for WNE and WC/WO stars including a 
spherically expanding envelope. The tracks are shown referring 
to different effective temperatures (cf. Sect. 2.1 for definitions) 
corresponding to the radius i?^, and the subphotospheric radii 
Rt=20 and Rt=io (same symbols as in Fig. 11). The shaded 
areas denote the WNE and WC/WO sequences respectively, 
which are obtained from different initial masses. The width of 
the sequences has been chosen to encompass the temperatures 
between T20 and Tio 

also of concern here. Again, in the WC/WO phases, the 
interior structure is not modified by the presence of the 
outfiowing envelope. 

The HR-diagram covering the WNE and WC/WO 
phases of the 85 Mq track is plotted in Figure 12. Shown 
are the different effective temperatures corresponding to 
the radius and the subphotospheric radii i?r=20 and 
Rj-io plotted in Fig. 11. Shaded areas indicate the posi- 
tion of the WNE and WC/WO sequences which are ob- 
tained for other initial masses. The width of these areas 
indicate the domain between T20 and Tiq. 

As a comparison with previous work we also plot the 
HR-diagram referring to T2/3 (Fig. 13). Accounting only 
for electron scattering yields the dashed line. The solid 
line shows the track obtained according to the treatment 
of Schaller et al. (1992) including electron scattering and 
line opacity in the wind, but using the observed wind ve- 
locities from Sect. 2. If instead we use a constant value of 
t^oo = 2200 km s~^ as in the original treatment of Schaller 
et al. we obtain the dotted line. Despite this difference the 
predictions for 72/3, and in particular the slopes in the HR 
diagram agree very well. This can easily be explained by 
the relative importance of electron scattering and the lines 
to the total opacity, as well as their different dependence 
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Fig. 13. HR-diagram for WNE and WC/WO stars including a 
spherically expanding envelope. Temperatures refer to an op- 
tical depth of r = 2/3. Solid line: preferred model adopting 
the Schaller et al. (1992) treatment including electron scatter- 
ing and line opacity, and observed wind velocities (see Sect. 2). 
Dotted line: same as solid line but using Voo = 2200 km 
following Schaller et al. Dashed line: T2/3 obtained using ob- 
served values for Voo but including only electron scattering. The 
shaded areas show the WNE and WC/WO sequence derived 
from the preferred model 

on M and t;oo (see Schaller et al. 1992). The shaded ar- 
eas in Fig. 13 denote the WNE and WC/WO sequence as 
obtained from our preferred treatment (solid line). 



6.1. Importance of CO opacities 

In this context it is worth mentioning briefiy the role of 
appropriate opacities for WC/WO stars. All evolutionary 
calculations for WR stars based on the recent OPAL data 
use opacities calculated with relative metal abundances 
corresponding to solar values. For evolved phases, in par- 
ticular for WC stars, which show the products of He burn- 
ing on their surface, this is clearly not valid anymore. We 
have therefore used the recent OPAL opacities for CO 
rich mixtures (Iglesias & Rogers 1993), which allow to 
correctly take the modified abundances of C and O into 
account. For typical WC model this results in a modest 
opacity enhancement (^ 20 %) in the outer layers, which 
yields a radius increase of up to 10 % only for (see 
Fig. 11). The He-burning lifetime of the 85 Mq model 
is, however, increased by only 0.32 %. For most purposes 
standard opacities are therefore clearly sufficient. 



Basic to our treatment is the assumption of radiative equi- 
librium in the expanding envelope. Although the same as- 
sumption is also made in non-LTE atmosphere models, 
which are used for comparisons in this work, one should 
be aware of a potential difficulty. Considering only the 
mechanical energy flux (neglecting the potential energy), 
we note that for the WN stars analysed by Hamann et 
al. (1993) this makes up between 0.4 % to 60 % of the 
radiative luminosity for the most extreme case. The influ- 
ence of this assumption on the derived stellar parameters 
is presently unknown. With the revision of luminosities 
(upward) and mass loss rates (downward) suggested by 
Schmutz (1995) the contribution of the mechanical flux to 
the total luminosity is, however, drasticly reduced. 

For the most luminous WNE and WC models from 
our evolutionary tracks, the kinetic energy flux may also 
reach up to ~ 20 % and 8 % of the radiative flux respec- 
tively. Accounting for the mechanical energy flux and also 
work done against the gravitational potential the lumi- 
nosity should be decreased by an amount raging between 
~ 0.14 (0.25) dex at logL/LQ = 6. and 0.02 (0.08) dex 
at logL/LQ = 5.5 referring to optical depths of r = 10 
(r = 2/3) respectively. As expected, differences are only of 
some importance at the highest luminosities, correspond- 
ing to M ^ lO-*M0yr-i. 

7. Comparison with observations 

A detailed comparison between predictions from evolu- 
tionary models and observations of WR stars, including 
luminosities, wind corrected temperatures (12/3), abun- 
dances, and population statistics has recently been pre- 
sented by Maeder & Meynet (1994). As shown above, the 
inclusion of the expanding envelope yields, for all cases 
except maybe WNL stars (cf. Sect. 4), the same results as 
the models used by Maeder & Meynet. Thus their con- 
clusions also apply to our models. On the other hand ba- 
sically two new predictions are made in this work. These 
will briefly be compared to observations in this paragraph. 

7.1. WN stars 

In Sect. 4 we have suggested the existence of a region lo- 
cated at 4.4 <, logT^ ^ 4.65, where WNL models are 
temporarily concentrated before evolving to higher tem- 
peratures. Interestingly the strip suggested by the mod- 
els with strong backwarming is found to be at the same 
"core" temperatures, where most observed WN stars still 
having hydrogen are placed by non-LTE spectral analysis 
(Hamann et al. 1993, Crowther et al. 1995b). Whether this 
is a pure coincidence is not clear. Statistical arguments 
about the relative WN population at different tempera- 
tures are surely not yet possible. On the other hand the- 
oretically there is no relation between the value of the H 
surface abundance and the proposed "WNL strip" , which 



show hydrogen on their surface (cf. Fig. 9). The location 
of the strip is only determined by the temperature. It is 
however clear, that the motion towards larger tempera- 
tures is related with an increasing (decreasing) He (H) 
abundance. Therefore stars located in the strip, and even 
more redwards of it, always have a non-zero H abundance. 
This seems to be in disagreement with some WNE-w stars 
analysed by Hamann et al. (1993). 

Regarding the extension of WNL stars, except 
WN-|-abs subtypes, we note that Crowther et al. (1995b) 
flnd values of R'j/s/R^ ~ 1.2-2.1 in contrast with 
R'j/s/R-k ^ 1-1 from Hamann et al. (1993), who claim 
that the continuum emerges from nearly hydrostatic lay- 
ers. Since both authors, however, use the same wind veloc- 
ity law with /3 = 1, this indicates that even for WNL stars 
the deflnition of the core radius depends on the adopted 
structure in the low velocity region, which is hardly con- 
strained by the observations. For WN stars with strong 
lines the ambiguous assignment of core radii from non- 
LTE models was already shown by Hillier (1991; see also 
Schmutz et al. 1992). This must be taken into account 
when comparing core radii and temperatures. 

7.2. WC stars 

As a second point our new calculations predict, for the 
flrst time, what could be called sub-photospheric radii, 
i.e. wavelength independent radii for WR stars with strong 
mass outflows. These radii, as well as the associated effec- 
tive temperatures, should allow a comparison to the so- 
called core radii determined from non-LTE spectral anal- 
yses of observed WC stars. As mentioned above, it has, 
however, to be kept in mind that even for non-LTE anal- 
yses the deflnition of a core radius may not be unique 
(Hillier 1991, Schmutz et al. 1992). 

We compare our results with the analysis of 25 Galactic 
WC stars of Koesterke & Hamann (1994). In their method 
two parameters are determined from line flts, taken apart 
the abundances. The basic parameters are the "trans- 
formed radius" 
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which is a measure of the wind density, and the effec- 
tive temperature T^. The corresponding "core" radius 
refers to depths of rRoss~ 10-20. To derive the other phys- 
ical quantities (R^, L^, M) requires the knowledge of the 
absolute visual magnitude, and the terminal velocity. The 
most direct comparison with their results is thus made in 
the T-f^-Rt plane. 

The evolutionary track of the 85 Mq model in the T^- 
Rt plane is shown in Fig. 14. We have plotted the track re- 
ferring to both the subphotospheric radii i?io and R20, but 
also using the radius of the wind-free hydrostatic model. 
For a comparison we have added a track where the wind 
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Fig. 14. Evolutionary track for WC/WO stars in the T-Rt 
plane. As in Figs. 11 and ?? the effective temperatures and 
radii refer to r = 10 (short-dashed), 20 (dotted), and the 
wind-free hydrostatic interior (solid). The dashed-dotted line 
is obtained with a wind opacity of 2(7e for r = 10. The WC-w 
and WC-s groups from the non-LTE analyses of Koesterke & 
Hamann (1994) are indicated by the shaded areas. Regions 
where He is completely ionised or recombines to He'*' are sep- 
arated by the thin dotted line (cf. Schmutz et al. 1992) 

opacity is 2(Te. The evolution proceeds essentially from the 
left to the right (increasing T). Note that independently 
of their initial mass and metallicity all WC models pop- 
ulate the same line in T-Rt plane, since these quantities 
are mainly determined by the mass loss rate, and R^, 
which are all related to the stellar mass^ (cf. Schaerer & 
Maeder 1992). Depending on their initial mass and metal- 
licity the starting position of the WC phase on this line, 
and in particular the surface abundances, do however vary 
considerably. 

Shaded areas indicate the domains of the analysed WC 
5 to WC 7 stars from Koesterke & Hamann (1994), for 
which they obtain good fits. Obviously the WC-w stars 
showing low wind densities and cooler effective tempera- 
tures are not matched by our models. If their continua 
should indeed arise from almost hydrostatic layers, as 
claimed by Koesterke & Hamann, the difference with our 
results cannot be due to an underestimate of the wind 
opacity, which would otherwise lead to lower temperatures 
and larger radii. 

On the other hand we note that our predictions refer- 
ring to r = 10 partly cover the same region as the WC 
5 to 7 subtypes with strong lines (WC-s). Subtypes ear- 

The obvious assumption is that the wind parameter /? and 
the opacity k are unchanged. 



higher temperatures. With ongoing evolution (C-|-0)/He 
increases, which implies that on one particular track later 
WC subtypes (cf. Smith k Maeder 1991) should be lo- 
cated at lower temperatures. Surprisingly Koesterke & 
Hamann indicate the opposite trend for WC 5 and 6 sub- 
types. As explained above, models with various metallic- 
ities and initial masses show, however, important surface 
abundance variations in the T-Rt diagram. Therefore one 
theoretically expects to find different WC subtypes mixed 
up at the same position in this diagram. 

In short we find that our evolutionary models, tak- 
ing into account the strong mass outfiow, roughly repro- 
duces the wind densities and effective temperatures of ob- 
served WC-s stars. Several points, however, show still dif- 
ferences with the results from non-LTE spectral analyses. 
In particular the existence of WC stars with "weak winds" 
(WC-w) at relatively low temperatures (T^ ~ 50 kK) can- 
not be understood with present day evolutionary mod- 
els. Particularly striking is also the large spread of carbon 
abundances derived by Koesterke & Hamann (1994) for 
a given WC subtype, contrarily to several previous deter- 
minations (cf. Smith & Hummer 1988, Eenens & Williams 
1992, Smith & Maeder 1991), which have also led to a sim- 
ple explanation of the observed subtype gradient (Smith 
& Maeder). Whether this may be related to the above- 
mentioned difficulties in comparing non-LTE model anal- 
ysis with evolutionary model predictions is, however, not 
clear yet. 

8. Conclusions 

A simple analytic method to describe the structure of a 
spherically expanding envelope with strong mass outfiow 
has been presented in this work. We solve for radiative 
equilibrium in the grey case and for mass conservation, 
while the momentum equation is replaced by adopting a 
parametrised velocity law. The structure is consistently 
connected to the hydrostatic stellar interior and allows us 
to study the infiuence of modified boundary conditions 
("backwarming") on the stellar structure and evolution 
(Sect. 2). On the same ground our treatment also permits 
to determine subphotospheric "core" radii and associated 
effective temperatures, which characterise the entire emer- 
gent spectrum of WR stars. 

We have applied the new method to evolutionary mod- 
els of WR stars. For WR stars showing still hydrogen on 
their surface (WNL's), models with strong backwarming 
predict a transient concentration in a strip located at tem- 
peratures 25000 ^ ^ 45000 K, indepently of the lumi- 
nosity. This is due to a small convective zone close to the 
surface, which is infiuenced by the outer boundary con- 
ditions and maintains a low density gradient, thereby al- 
lowing to temporarily stabilise the stellar radius on the 
blueward track (see Sect. 4). Surprisingly this strip coin- 
cides quite well with the "degenerate " location, where 



hand, if a moderate backwarming is adopted, the tracks 
in the HR-diagram are basically unchanged with respect 
to the usual treatment using plane parallel atmospheres, 
and the "WNL strip" disappears. 

In any case, including all the WR phases up to the end 
of He-burning, the interior structure and hence the evo- 
lution is never modified by the presence of an outfiowing 
envelope (see Sects. 4 to 6). It is well known that the dom- 
inant process determining the evolution in this domain is 
mass loss (cf. Chiosi & Maeder 1986). We have now shown, 
that it is indeed sufficient to consider only the effect of a 
progressive reduction of the stellar mass. The other conse- 
quences of mass loss, including backwarming from a spher- 
ically expanding envelope, have no effect on the deep inte- 
rior structure. This does, however, not generally hold for 
the stellar parameters (radii, effective temperatures) and 
in particular does not apply to the spectral appearance, 
which obviously strongly depends on the "details" of the 
mass loss phenomenon. 

For hydrogen free (WNE, WC and WO) stars with 
strong mass loss we have calculated subphotospheric radii 
corresponding to optical depths of r ~ 10-20 (Sects. 5 and 
6). The radius increase with respect to wind-free stars 
(both hydrostatic or hydrodynamic models) may be up 
to a factor of ~ 4 for the most luminous WNE or WC 
stars. Provided the hydrodynamic structures (particularly 
at large depths) adopted in both our approach and in 
non-LTE models do not differ significantly, our subphoto- 
spheric radii should be comparable with "core" radii de- 
rived from non-LTE spectral analyses. Comparisons with 
observations show a better agreement than previous mod- 
els (Sect. 7). Future progress awaits full hydrodynamic 
calculations of coupled interior and atmosphere models. 

Subphotospheric radii also allow a direct assignment 
of spectra to evolutionary tracks. Compared with wind- 
free evolutionary models this implies e.g. H (He) ionising 
fiuxes (derived from the non-LTE models of Schmutz et 
al. 1992), which are reduced by up to ~ 30 (60) %. Fi- 
nally our treatment of the envelope now permits to pursue 
the study of the spectral and interior evolution of mas- 
sive stars presented by Schaerer (1995) and Schaerer et 
al. (1995ab) to post-MS phases including WN, WC and 
WO stars. 

Last, but not least, we have pointed out the possible 
importance of the iron opacity peak for the acceleration of 
WR winds in the optically thick part (Sect. 4). This may 
well be essential for the understanding of the dynamics of 
Wolf-Rayet winds. 
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